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Executive  Summary 
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Mach  3.0  jet  were  compiled  and  submitted  for  publication.  Efforts  now  turn 
to  completing  laboratory  experiments  on  the  heated,  conic-section  nozzle, 
supersonic  jet  in  the  U.  Miss,  facility.  Initial  shadowgraph  images  provide 
insight  into  the  generation  of  shock-like  events  in  the  jet’s  near-field. 
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1  Project  Objectives  and  Status 

1.1  Review  of  Program  Objectives 

Broadly,  the  objectives  of  this  effort  focus  on  the  finite  amplitude  disturbances  generated 
by  supersonic  jets,  their  impact  on  measurement  and  modeling  of  the  jet  noise,  and  the 
importance  of  these  events  in  evaluation  of  any  jet  noise  control  strategy.  Experimental 
measurements  have  been  obtained  for  shock-free  and  shock-containing  supersonic  jets. 

•  An  acoustic  survey  of  the  sound  field  indicates  that  the  sound  source  is  in  the  post¬ 
potential  core  region.  Contours  of  constant  amplitude  are  curved  indicating  that  a 
distribution  of  convection  velocities  is  associated  with  the  source  disturbances. 

•  Evaluation  of  the  “shock”  content  in  the  pressure  time-series  indicates  finite  amplitude 
events  can  be  observed  immediately  in  the  near-held  as  well  as  being  the  result  of 
steepening  in  the  mid-held. 

•  Evaluation  of  the  decay  along  the  peak  emission  angle  indicates  a  transition  from 
cylindrical  to  spherical  spreading  occurs  between  12  and  33  jet  diameters  in  radius. 

•  Evaluation  of  the  time  series  of  the  pressure  along  the  peak  emission  angle  indicates 
some  waves  may  coalesce  while  others  may  steepen. 

•  Shadowgraph  images  verify  hnite  amplitude  disturbances  are  emitted  both  from  the 
near-nozzle  how  and  the  post-potential  core  region. 
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Figure  1.  The  acoustic  disturbances  in  the  near-field  of  the  jet  appear  as  conical-shaped, 
expanding  surfaces.  Flow  is  left-to-right. 


1.2  Project  Status 

The  results  of  the  most  recent  set  of  measurements  conducted  in  the  UT  Austin  facility 
have  been  submitted  to  the  AIAA  Journal  for  publication.  Now,  the  focus  of  the  remaining 
experimental  efforts  shift  to  the  U.  of  Miss,  facility  where  acoustic  and  flow  measurements 
will  be  conducted  with  a  hot-jet  whose  acoustic  Mach  number  is  higher  than  that  of  the  cold 
jet  used  for  the  UT  Austin  measurements. 

2  Activity  for  Current  Reporting  Period 

2.1  Disturbances  in  the  Near-Field  of  the  Jet 

Isosurfaces  of  the  density  were  analyzed  from  the  large-eddy  simulation  (LES)  of  the  over- 
expanded,  conic-section  nozzle  {NPR  =  3.96,  TTR  =  3.35,  Mj  =  1.55,  Ma  =  2.34).  The 
iicar-ficld  structure  is  characterized  by  conical  shaped  surfaces  as  illustrated  in  Figure  1.  The 
cone  half-angle  of  these  structures  is  associated  with  the  convection  velocity  of  the  underlying 
disturbance.  Convection  speeds  between  0. 7-0.8  of  the  jet’s  acoustic  Mach  number  suggest 
angles  between  32  and  37  degrees.  However,  this  velocity  is  not  constant  along  the  jet  axis. 

A  contour  plot  of  the  dilatation  from  the  LES  is  shown  in  Figure  2.  The  annotations 
(black  lines)  highlight  the  angle  that  the  disturbance  appears  to  make.  These  angles  are 
steeper  at  the  jet  exit  and  flatten  out  in  the  downstream  direction.  This  indicates  that  the 
underlying  disturbances  must  be  slowing  down  as  they  move  along  the  jet.  Certainly,  the 
far-field  directivity  has  a  particular  peak,  but  in  the  near-field  the  underlying  disturbances 
should  be  considered  to  result  from  some  distribution  of  convection  speeds. 

Shadowgraph  images  were  recorded  in  the  near-field  of  both  the  nozzle  exit  and  post¬ 
potential  core  regions  of  the  jet.  A  retro-reflective  shadowgraph  setup  was  used  to  record 
the  images.  The  near-nozzle  region  exhibits  a  large  number  of  shock-like  disturbances  being 
generated.  While  not  regular,  the  images  recorded  in  the  post-potential  core  region  also 
exhibit  a  stochastic  occurrence  of  finite  amplitude  disturbances.  The  shape  and  direction  of 
these  disturbances  suggest  that  they  originated  in  the  post-potential  core  region  of  the  jet. 
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(a)  Nozzle-Exit  Region  (b)  Post-Potential  Core  Region 


Figure  3.  Shadowgraph  imagery  showing  finite  amplitude  disturbances  emanating  from  both 
the  near- nozzle  region  and  the  post-potential  core  region.  The  vertical  line  in  (b)  is  located  8 
jet  diameters  downstream  from  the  nozzle  exit. 


This  observation  suggests  that  there  is  a  distribution  of  disturbances  along  the  jet  that  reach 
supersonic  phase  speeds  with  sufficient  amplitude  to  generate  finite  amplitude  waves  in  the 
near-field  of  the  jet. 

2.2  Wave  Coalescence 

An  inspection  of  Figure  3(a)  indicates  that  the  shock-like  structures  are  coalescing  as  they 
propagate  away  from  the  jet.  This  supports  concepts  suggested  by  the  acoustic  measurements 
performed  by  the  UT  Austin  group  (see  Appendix  A)  This  process  is  only  well  understood 
in  a  1-D  sense.  The  3-dimensional  merging  of  finite  amplitude  disturbances  is  not  well 
understood,  and  the  effect,  if  any,  on  the  jet’s  near-  end  far-field  directivity  is  unknown  at 
this  time. 

3  Technical/Cost  Status  &  Problem  Areas 

None. 

4  Publications,  Meetings,  and/or  Travel 

•  Dr.  Tinney’s  student,  Remain  Fievet,  presented  recent  experimental  results  at  the  2014 
Texas  Fluid  Dynamics  Conference.  The  abstract  is  included  in  Appendix  B. 

•  Results  from  recent  work  performed  by  Dr.  Tir.ney  and  his  students  have  been  sub¬ 
mitted  to  the  AIAA  Journal.  The  submission  draft  is  included  in  Appendix  A. 
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•  Dr.  Tinney’s  former  student,  Dr.  Baars,  completed  a  journal  submission  to  Physics 
of  Fluids  on  the  ‘TOD-Based  Spectral  Higher-Order  Stochastic  Estimation.”  This 
concept  was  included  in  the  original  proposal  for  this  project,  so  it  is  mentioned  here 
and  included  in  the  publications  list. 

4.1  Running  List  of  Publications  Produced 

[1]  Baars,  W..J.,  Tinney,  C.E.,  Murray,  N.E.,  .Jansen,  B.J.  &  Panickar.  P.  (2011)  ‘The  Effect  of 
Heat  on  Turbulent  Mixing  Noise  in  Supersonic  .Jets,”  in  4^^  AIAA  Aerospace  Sciences  Meeting^  Paper 
2011-1029  (Orlando,  FL). 

[2]  Baars,  W.J.,  Tinney,  C.E.  Sz  Wochner,  M.S.  (2012)  '‘Nonlinear  Noise  Propagation  from  a  Fully 
Expanded  Mach  3  .Jet,”  in  50th  AIAA  Aerospace  Science  Meeting^  Paper  2012-1177  (Nashville,  TN). 

[3]  Murray,  N.E.,  Lyons,  G.W.,  Tinney,  C.E.,  Donald,  B.,  Baars,  W..J.,  Thurow,  B.S.,  Haynes, 
R.H.  Sz  Panickar,  P.  (2012)  “A  Laboratory  Framework  for  Synchronous  Near/Far-Field  Acoustics  and 
MHz  PIV  in  High-Temperature,  Shock-Containing.  .Jets,”  in  Proceedings  of  the  Intemoise  2012/ASME 
NCAD  Meeting,  ASME/NCAD-1270  [invited]  (New  York  City,  NY). 

[4]  Baars,  W..J.  &  Tinney,  C.E.  (2012)  “Scaling  Model  for  Nonlinear  Supersonic  Jet  Noise,”  in  Bidletin 
of  the  American  Physical  Society,  57:17,  Abstract  D24:8  (San  Diego,  CA). 

[5]  Panickar,  P.,  Erwin,  J.,  SInha,  N.,  Murray,  N.E.  k,  Lyons,  G.W.  (2013)  “Localization  of 
Acoustic  Sources  in  Shock-Containing  Jet  Flows  Using  Phased  Array  Measurements,”  in  5P^  AIAA 
Aerospace  Sciences  Meeting,  2013-0613  (Grapevine,  TX). 

[6]  Haynes,  R.H.,  Brock,  B.A.  Iz  Thurow,  B.S.  (2013)  "Application  of  MHz  Frame  Rate,  High  Dy¬ 
namic  Range  PIV  to  a  High-Temperature,  Shock-Containing  Jet,”  in  5P^  AIAA  Aerospace  Sciences 
Meeting,  2013-0774  (Grapevine,  TX). 

[7]  Fievet,  R.,  Tinney,  C.E.,  Murray,  N.E.,  Lyons,  G.W.  ^  Panickar,  P.  (2013)  “Acoustic  Source 
Indicators  using  LES  in  a  Fully  Expanded  and  Heated  Supersonic  Jet,”  in  19^^^  AIAA/CEAS  Aeroa- 
coustics  Conference,  2013-2193  (Berlin,  Germany). 

[8]  Baars,  W.J.  Sz  Tinney,  C.E.  (2013)  “Quantifying  Crackle-Inducing  Acoustic  Shock-Structures  Emit¬ 
ted  by  a  Fully- Expanded  Mach  3  Jet,”  in  19^^  AIAA/CEAS  Aeroacoustics  Conference,  2013-2081 
(Berlin,  Germany). 

[9]  Baars,  W.J.  Sz  Tinney,  C.E.  (2013)  “A  Temporal  and  Spectral  Quantification  of  the  Crackle  Com¬ 
ponent  in  Supersonic  Jet  Noise,”  in  Symposium  on  Fluid- Structure- Sound  Interactions  and  Control. 
p.  170  (Shatin,  Hong  Kong). 

[10]  Baars,  W.J.,  Tinney,  C.E.  Sz  Wochner,  M.S.  “Nonlinear  Distortion  of  Acoustic  Waveforms  from 
High-Speed  Jets,”  J.  of  Fluid  Mechanics  [in  review]. 

[11]  Baars,  W.J.  (2013)  Acoustics  from  High-speed  Jets  with  Crackle,  Ph.D.  Dissertation,  The  University 
of  Texas  at  Austin. 

[12]  Fievet,  R-.,  Baars,  W.J.,  Silva,  D.  Sz  Tinney,  C.E.  (2013)  “High  Fidelity  Measurements  in  the  Far- 
Field  of  a  Mach  3  Jet,”  in  Bulletin  of  the  American  Physical  Society,  58:18,  Abstract  D24:7  (Pittsburgh, 
PA). 

[13]  Brock,  B.A.,  Haynes,  R.H.,  Thurow,  B.S.,  Lyons,  G.W.  Sz  Murray,  N.E.  (2014)  “An  Exam¬ 
ination  of  MHz  Rate  PIV  in  a  Heated  Supersonic  Jet,”  in  AIAA  Science  and  Technology  Forum  and 
Exposition  (SciTech  2014)^  2014-1102  (National  Harbor,  MD). 

[14]  Baars,  W.J.  Sz  Tinney,  C.E.  (2014)  “Temporal  and  Spectral  Quantification  of  the  Crackle  Compo¬ 
nent  in  Supersonic  Jet  Noise,”  in  Y.  Zhou,  Y.  Liu,  L.  Huang  Sz  D.H.  Hodges,  eds.,  Fluid- Structure- Sound 
Interactions  and  Contjvl,  Lecture  Notes  in  Mechanical  Engineering,  pp.  205-210  (Springer  Berlin  Hei¬ 
delberg)  . 
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[15]  Baars,  W.J.  Sz  Tinney,  C.E.  (2014)  ‘  POD-Based  Spectral  Higher-Order  Stochastic  Estimation,” 
Physics  of  Fluids  [in  review]. 

[IG]  Fievet,  R.  &  Tinney,  C.E.  (2014)  ‘‘Nonlinear  Acoustic  Distortion  in  Supersonic  Lab-Scaled  Jet 
Vicinity,”  in  Texas  Fluid  Dynamics  Meeting. 

[17]  Fievet,  R.,  Tinney,  C.E..  Baars.  W.J.  &  Hamilton,  M.  (2014)  “Acoustic  Waveforms  Produced 
by  a  Laboratory  Scale  Supersonic  Jet,”  AIAA  J.  [in  review). 

5  Planned  Activities  for  Next  Reporting  Period 

The  remaining  laboratory  activities  at  U.  Miss,  will  be  complete  during  the  next  reporting 
period.  These  measurements  will  include  a  survey  of  the  sound  field,  line-array  microphone 
measurements  of  the  decay  in  the  peak  noise  direction,  LDA  measurements  of  the  nozzle  exit 
turbulent  boundary  layer  (relevant  for  analysis  of  instability  waves  of  the  jet),  and  additional 
shadowgraph  measurements. 
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Appendix  A:  Submitted  to  AIAA  Journal  for  publication: 

Acoustic  Waveforms  Produced  by  a  Laboratory  Scale  Supersonic  Jet 
Fibvc't,  Tiinicy,  Baars  &  Hamilton 
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Acoustic  waveforms  produced  by  a  laboratory  scale 

supersonic  jet. 

Romain  Ficvct  and  Charles  E.  Tinney 

Center  for  Aeromechanics  Research,  The  University  of  Texas  at  Austin,  Austin,  TX  18112,  USA 

Woutijri  J.  Baars 

The  University  of  Melbourne,  Parkville,  Victoria  3010,  Australia 
Mark  F.  Hamilton 

Applied  Research  Laboratories,  The  University  of  Texas  at  Austin,  Austin,  TX  18113,  USA 

The  spatial  evolution  of  acoustic  waveforms  produced  by  a  Mach  3  jet  are  investi¬ 
gated  using  both  1/4  inch  and  1/8  inch  pressure  field  microphones  located  along  rays 
emanating  from  the  post  potential  core  where  the  peak  sound  emission  is  found  to 
occur.  The  measurements  are  acquired  in  a  fully  anechoic  chamber  where  ground,  or 
other  large  surface  reflections  are  minimal.  The  calculation  of  the  OASPL  along  an  arc 
located  at  95  jet  diameters  using  120  planar  grid  measurements  are  shown  to  collapse 
remarkably  well  when  the  arc  array  is  centered  on  x/Dj  =  17.5  which  is  where  the  most 
intense  sound  is  generated.  Various  statistical  metrics,  including  the  quadrature  spec¬ 
tral  density,  number  of  zero  crossings,  the  skewness  of  the  pressure  time  derivative  and 
the  integral  of  the  negative  part  of  the  quadrature  spectral  density,  are  exercised  along 
the  peak  emission  path.  These  metrics  are  shown  to  undergo  rapid  changes  within  2 
meters  from  the  source  regions  of  this  laboratory  scale  jet.  The  sensitivity  of  these  find¬ 
ings  to  both  transducer  size  and  humidity  effects  are  discussed.  A  visual  extrapolation 
of  these  nonlinear  metrics  toward  the  jet  shear  layer  suggests  that  these  waveforms 
are  initially  skewed  at  the  source.  The  agreement  between  a  wave  packet  modeling 
of  the  jet  near-field  and  the  experiment  allows  identifying  the  transition  zone  from 
cylindrical  to  spherical  pressure  decay  law,  and  evaluate  the  source  intensity,  which  is 
used  to  estimate  the  shock  formation  distance  and  the  Gol’dberg  number. 
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A 

Nomenclature 

=  non-dimensional  wave  packet  parameter,  L/Dj 

aj 

—  speed  of  sound  of  the  jet  at  the  exit  of  the  nozzle,  in/s 

O-cjo 

—  ambient  speed  of  sound,  m/s 

Dj 

=-  jet  exit  diameter,  m 

Ic 

=  cut-off  frequency,  Hz 

fj 

=  jet  characteristic  frequency,  Hz 

fo 

=  peak  spectral  frequency  along  the  propagation  path,  Hz 

Gpp 

=  one-sided  power  spectral  density  of  the  pressure,  Pa^/Hz 

L 

=  wave  packet  model  scale 

ic 

—  absorption  length  scale 

Ma 

=  acoustic  Mach  number,  Uj/aao 

Mi 

—  gas  dynamic  Mach  number  of  the  jet  at  the  exit  of  the  nozzle,  f/j/oj 

Me 

=  Convective  Mach  number,  Vd^oo 

N 

—  Number  of  microphones 

N,c 

=  Number  of  zero  crossings,  #/s 

V 

=  pressure  waveform,  Pa 

P 

=  time  derivative  of  the  pressure  waveform.  Pa  s 

Po 

=  total  pressure  of  the  jet  at  the  exit  of  the  nozzle,  K 

Poo 

=  ambient  pressure,  K 

Q 

=  normalized  quadrature  spectral  density,  1  /Hz 

Qneg 

=  integral  of  the  negative  part  of  Q 

Qp^p 

—  quadrature  spectral  density  of  the  pressure,  Pa^  /Hz 

T 

—  radial  coordinate  relative  to  the  jet  centerline. 

Tc 

=  cylindrical  shock  formation  distance,  in 

Ts 

=  spherical  shock  formation  distance,  m 

%RH 

=  percent  relative  humidity 

1 

—  skewness  factor 

Spp 

=  two-sided  power  spectral  density  of  the  pressure,  Pa^/Hz 

Stoj 

—  Strouhal  number:  JDj/Uj 

-V 

=  standard  deviation  of  p(t) 

t 

—  time,  s 

To 

=  total  temperature  of  the  jet  at  the  exit  of  the  nozzle,  K 

Tj 

—  static  temperature  of  the  jet  at  the  exit  of  the  nozzle,  K 

=  ambient  temperature,  K 

T 

=  retarded  time,  s 

Uc 

—  convective  speed  of  the  jet,  m/s 

Uj 

=  velocity  of  the  jet  at  the  exit  of  the  nozzle,  m/s 

w 

=  wave  packet  amplitude 

w 

' ' error 

—  Wave  packet  error 

X 

=  axial  coordinate  relative  to  the  jet  exit  plane 

X 

—  plane  wave  shock  formation  distance 

a 

—  absorption  factor,  Np/m 

=  coefficient  of  nonlinearity 

r 

=  Gol’dberg  number  for  plane  waves 

Ac 

=  modified  cylindrical  Gol’dberg  number 

A, 

■=  modified  spherical  GoPdberg  number 

<!> 

—  Mach  wave  angle  measured  from  the  jet  centerline,  deg. 

p 

—  distance  along  a  ray  emanating  from  the  peak  emission  location  of  the  jet,  m 

Poo 

=  ambient  density,  kg/m^ 

^p,0 

—  pressure  source  strength,  Pa 

e 

—  angle  measured  from  the  jet  centerline,  deg. 

^pp 

—  pressure  cross  correlation,  Pa^ 

FSS 

=  Fine  Scale  Similarity  spectra 

LSS 

—  Large  Scale  Similarity  spectra 

NPR 

=  Nozzle  pressure  ratio,  Pol  Poo 

OASPL 

f  —  Overall  sound  pressure  level,  dB,  ref;20/xPa^/Hz 

PDF 

—  Probability  density  function. 

4^  AiAA  1 - - 2  rs - r> - *  .4.^ 
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I.  Introduction 

When  the  eonveetivc  speed  of  large  seale  turbulent  structures  travel  faster  than  the  sound  speed 
of  the  surrounding  gas,  they  emit  Maeh  waves  which  propagate  at  angles  relative  to  the  jet  axis 
determined  by  =  eos^^{aoo/Uc)  |1  6|.  Because  these  Maeh  waves  arc  generated  by  turbulence, 
they  arc  classified  as  a  component  of  turbulent  mixing  noise,  which  is  the  most  significant  source  of 
noise  for  jets  operating  at  supersonic  acoustic  Mach  numbers.  If  the  amplitude  of  the  Mach  wave  is 
large  enough  to  overcome  viscous  relaxation  effects  |7,  8],  the  waveform  will  steepen  with  distance 
from  its  source,  which  is  more  easily  observed  imdcr  full-scale  conditions  where  the  characteristic 
frequencies,  and  hence  viscous  relaxation  effects,  arc  low  |9|, 

While  sawtooth-like  waveforms  have  been  observed  in  both  laboratory  and  full  scale  jet  noise 
studies,  the  source  and  nature  of  these  shapes  remains  unknown.  Most  laboratory  and  full-scale 
studies  are  conducted  using  arc  arrays  of  microphones  placed  in  the  far-ficld  of  the  jet  (beyond  60 
or  so  jet  diameters).  And  so  the  shape  and  evolution  of  these  Maeh  waves,  as  they  propagate  from 
the  source  regions  of  the  flow  to  the  far-field,  are  not  well  understood,  A  recent  study  by  Baars, 
Tinncy  Woehner  (2012)  |10|  mapped  the  acoustic  field  of  an  unheated  supersonic  jet  (between 
x/Dj  =  [5  :  10  :  145)  and  r/Dj  =  [25  :  10  :  95])  operating  at  a  gas  dynamic  Maeh  number  of  3 
to  unveil  the  sound  field  produced  by  jets  comprising  supersonic  acoustic  Mach  numbers.  Baars  & 
Tinney  (2014)  |6|  later  postulated  that  steepened  waveforms  are  produced  either  at  the  source  or 
develop  rapidly  in  the  ncar-ficld  through  cumulative  nonlinear  distortion  effects  (much  earlier  than 
where  typical  far-ficld  measurements  are  acquired).  Recent  Large- Eddy-Simulations  by  Nichols  et  aL 
(2013)  |11|  showed  sawtooth  waves  being  emitted  from  within  the  jet  shear  layer  and  that,  despite 
further  steepening  from  cumulative  nonlhicar  distortions,  it  was  likely  that  steepened  waveforms 
were  also  being  generated  within  the  jet,  Anderson  Sc  Frcimd  (2012)  |12|  observed  similar  features 
using  direct  numerical  simulations  and  went  on  to  illustrate  the  clumping  together  of  shorter  waves 
in  the  very  near  field  of  supersonic  mixing  layers.  Such  observations  would  suggest  that  large 
amplitude  waves  do  not  form  on  their  own,  but  rather  form  by  partnering  with  other  smaller  waves 
to  produce  larger  and  more  steepened  waves  further  from  the  jet.  Thus,  the  observations  of  Nichols 
et  aL  (2013)  |11|,  Baars  Sc  Tinney  (2014)  |6|  and  Anderson  Sc  Freund  (2012)  [12]  suggest  that  the 
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cinulsion  of  skewed  waveforms,  from  within  the  jet,  proeeeded  by  either  rapid  waveform  steepening 
or  the  partnering  of  similar  waves  along  the  propagation  path  are  responsible  for  sawtooth  like 
struetures  observed  in  the  far-held.  And  so,  an  important  question  that  we  seek  to  address  here 
is  whether  the  aeoustie  waveforms,  produeed  by  Maeli  waves  emanating  from  laboratory  seale  jets, 
undergo  eumulative  nonlinear  distortions,  or  steepen  through  other  meehanisms. 

Here,  we  will  replieate  portions  of  the  experiment  performed  by  Baars  et  ai  (2012)  [10]  who 
used  1  /4  ineh  inierophones  to  eharaeterize  the  aeoustie  field  of  a  Mach  3  jet.  Contrary  to  that  study, 
we  will  employ  higher  fidelity  instruments  (1/8  inch  pressure  field  mierophones),  positioned  closer 
to  the  jet  flow  and  along  lines  emanating  from  the  post  potential  core  region  where  the  prominent 
source  of  sound  is  known  to  reside.  As  was  done  before,  various  well  known  statistical  metrics 
for  quantifying  waveform  steepening  and  coalescenee  will  be  exercised  to  determine  the  shape  and 
nature  of  the  waveforms  as  they  propagate  along  rays  centered  on  the  post-potential  core.  The 
effect  of  absorption  on  these  various  metrics  will  be  studied,  as  the  measurements  of  Baars  et  aL 
(2012)  [10]  were  conducted  at  higher  relath^e  humidities  than  the  measurements  acquired  with  the 
1  8  inch  microphones.  Doing  so  will  increase  the  likelihood  of  capturing,  over  a  broader  spectral 
range,  distortions  in  the  acoustic  waves  as  they  radiate  to  the  far-field  of  laboratory  scale  supersonic 
jets. 

II.  Experimental  Description 

A.  Facility  overview 

Both  the  experiments  of  Baars  et  aL  (2012)  [10]  and  the  current  higher  fidelity  measurements 
were  conducted  in  the  fully  anechoic  chamber  and  open  jet  wind  tunnel  at  The  University  of  Texas 
at  Austin.  An  illustration  of  the  chamber  is  shown  in  figure  1  and  comprises  internal  dimensions 
(wedge-tip  to  wedge-tip)  of  5.74  m(lcngth)x4.52  m(widtli)x3.66  m(height).  The  acoustic  wedges 
that  arc  used  to  produce  this  echo  free  environment  are  composed  of  46  cm  of  fire  retardant  melamine 
foam  backed  by  a  10  cm  air  cavity  followed  by  14  cm  of  recycled  cotton  fibre.  The  findings  from 
impedance  tube  tests  performed  at  ETS  Lindgrcn  (Cedar  Park,  Texas)  arc  shown  in  figure  2  to 
demonstrate  how  this  particular  wedge  design  produces  a  normal  incidence  absorption  coeffieient  of 
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99%  above  100  Hz  for  an  air  cavity  depth  of  10  cm;  the  small  dip  around  250  Hz  is  due  to  the  flow 
resistive  properties  of  the  melamine  foam. 

A  jet  rig  is  positioned  along  the  eenterline  of  the  ehaiiiber  and  tunnel,  as  shown  in  figure  1  with 
a  1.49  m^  square  aperture  on  the  back  wall  to  allow  entrained  flow  to  enter  into  the  chamber.  Both 
the  jet  and  entrained  flows  exhaust  downstream  through  a  1 .83  m^  aeoustieally  treated  square  duct 
and  with  the  assistance  of  a  vane  axial  fan  to  reduce  recirculations  from  forming  inside  the  chamber. 
The  nozzle  tested  is  the  same  £is  that  described  by  Baars  et  al  (2012)  [10]  and  was  designed  to 
produce  a  perfectly  expanded  flow  at  a  gas  dynamic  Mach  number  of  3.0  with  an  exit  diameter 
of  Dj  =  0.0254  m  using  unheated  air.  Nozzle  pressure  ratios  (NPR)  are  monitored  and  controlled 
using  a  National  Instruments  ConipactRIO  system  which  is  capable  of  reducing  variations  in  the 
set-point  NPR  to  less  than  1%  over  extended  periods  of  time.  Additional  details  concerning  this 
facility  eaii  be  found  in  Donald  et  al,  (2014)  [13]  and  Baars  (2013)  (14]. 


Fig.  1:  The  fully  aiiechoic  chamber  at  UT-Austin  during  the  current  set  of  measurements. 

B.  Jet  operating  conditions 

Five  separate  tests  were  conducted  in  September  2013  over  a  duration  of  2  days  during  which 
atmospheric  conditions  were  carefully  recorded;  see  table  1.  While  the  fully  anechoic  chamber  and 
wind  tunnel  are  located  in  a  large  warehouse;  outside  air  is  drawn  in  to  prevent  the  building  from 
collapsing  during  runs.  Therefore,  the  environment  inside  the  chamber  is  that  of  the  air  outside  of  the 
building.  The  nozzle  was  operated  at  a  nozzle  pressure  ratio  of  36.7  which  was  expected  to  produce 
a  perfectly  expanded  flow  at  the  exit  of  the  nozzle  (based  on  quasi  1-D  isentropic  compressible 
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Fig.  2:  Normal  incidence  absorption  coefficient  of  the  UT  Austin  chamber  wedge  with  10  cm  air 

cavity. 

flow’  equations).  Shadowgraph  images  later  revealed  the  presence  of  weak  shock  structures  at  this 
pressure  ratio.  How^ever,  no  broadband  shock  tones  or  screech  were  observed  in  any  of  the  acoustic 
spectra,  and  so  these  shocks  were  not  expected  to  have  any  influence  on  the  measurements  reported 
here  (or  those  acquired  in  2011  with  the  1/4  inch  microphones  [6,  10,  14]).  A  common  set  of 
values  for  Poo,  T^oo  and  To  W’cre  chosen  by  averaging  over  the  test  conditions  reported  in  table  1. 
The  resultant  differential  error  bctw’ccn  these  tests  was  found  to  be  within  1.5%  of  the  average 
value.  Included  in  this  table  are  the  average  values  reported  by  Baars  et  al.  (2012)  [10],  (shown  in 
parenthesis)  which  compare  reasonably  well  to  the  conditions  reported  here.  The  largest  diserepaney 
between  the  current  set  of  measurements  and  the  measurements  of  Baars  et  ai  (2012)  [10]  reside  in 
the  atmospheric  conditions  (relative  humidity  and  atmospheric  temperature),  with  percent  relative 
humidities  of  75.5,  63.2  and  48.8  being  recorded  by  Baars  et  aL  (2012)  [10]  on  the  three  days  of 
testing,  respectively. 

The  convective  speed  of  the  sound  producing  events  is  an  important  quantity  and  is  assumed 
to  be  O.SUj  [4].  This  results  in  Mach  wave  angles  around  0o.8  =  45°,  which  coincide  quite  well  with 
the  peak  noise  emission  angle  shown  in  figure  3;  this  angle  is  determined  to  be  where  the  overall 
sound  pressure  level  (OASPL)  contour,  also  shown  in  figure  3  from  the  (rr,  r)  grid  nieasurcments  of 
Baars  et  ai  (2012)  [lO],  is  iiiaximum.  The  lobe-shape  pattern  in  the  OASPL  is  manifest,  owing 
to  the  eoinpetiiig  effects  of  eonveetioii  and  refraction  that  is  inherent  to  jet  noise.  However,  closer 
inspection  of  the  OASPL  ridge  demonstrates  how  contour  levels,  residing  closer  to  the  jet,  manifest 
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Tabic  1:  Testing  conditions  for  the  1/8  inch  microphone  study.  Average  values  from  the  1  '4  hicli 
study  arc  shown  in  parenthesis.  Note:  the  microphone  locations  corresponding  to  the  45°  and  455 

arrays  are  provided  in  table  2. 


Array  angle 

Measured  values 
35° 

45? 

455 

52° 

65° 

Average 

Poo,  kPa 

Po/Poo 

100.0 

100.1  100.0 

36.7  ±1%  (±4.5%) 

100.2 

100.2 

100.1  (100.9) 
36.7 

Too,K 

302 

304 

306 

303 

303 

304  (289) 

To,K 

301 

298 

301 

304 

304 

302  (288) 

%RH 

43 

Computed  values 

51 

40 

56 

58 

N'A 

Ti,K 

108 

106 

108 

109 

109 

108  (103) 

Gj,  m  s 

208 

207 

208 

209 

209 

208  (203) 

Uj,  m/s 

624 

620 

624 

627 

627 

624  (609) 

fj,  kHz 

24.6 

24.4 

24.6 

24.7 

24.7 

24.6  (24.0) 

Mj  =  Uj/aj 

3.0 

3.0 

3.0 

3.0 

3.0 

3.0 

a  ™  t/7  ! floo 

1.78 

1.78 

1.78 

1.80 

1.80 

1.79  (1.79) 

Me  =  O.Sf/j/a^o 

1.43 

1.42 

1.42 

1.44 

1.44 

1.43  (1.43) 

deg. 

45 

45 

45 

46 

46 

45  (46) 

deg. 

37 

36 

37 

37 

37 

37  (37) 

Mach  waves  with  angles  (0)  coinciding  with  convective  speeds  on  the  order  of  0.7Uj  (identified  in 
figure  3  by  0o.7  =  35°).  Contrarily,  in  regions  further  from  the  jet,  Mach  wave  angles  coincide  with 
structures  moving  at  O.SUj  (identified  in  figure  3  by  0o.8  =  45°).  Variations  in  the  Mach  wave  ridge 
are  attributed  to  differences  in  the  convective  speed  of  the  more  pronounced  sound  producing  events 
and  their  locations  in  the  post  potential  core  region  relative  to  one  another. 


C.  Instrumentation  and  microphone  placement 

Because  results  from  two  different  studies  employing  different  kinds  of  instruments  are  being 
evaluated,  the  details  concerning  both  are  provided  here  for  coniplctcncss.  The  first  of  the  acoustic 
measurements  is  described  by  Baars  et  al.  (2012)(10|  and  was  performed  in  December  of  2011  using 
four  1  4  inch  prcpolarizcd,  pressure- field  condenser  microphones  (PCB  model  377B10  capsules  with 
426B03  preamplifiers)  positioned  at  grazing  incidence  to  the  jet  axis.  These  lEPE  type  transducers 
were  powered  using  a  NI  PXI-4472  card  which  then  digitized  the  signals  at  102.4  kS/s  with  24 
bit  resolution  (ilOv)  for  a  minimum  of  2^^  samples.  Signal  attenuation  occurs  at  0.84  times  the 
Nyquist  frequency  with  the  built  in  anti-aliasing  filters  on  the  NI  PXI-4472  card.  The  placement  of 
these  microphones  involved  both  a  2-D  grid  of  points  in  the  (x,  r)  plane  as  well  as  a  line  array  along 
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the  peak  emission  path,  which  was  found  to  emanate  from  a  region  centered  around  x/Dj  =  17.5. 
Spectral  based  metrics  arc  constructed  using  2^^  samples  per  block,  thus  resulting  in  a  frequency 
resolution  of  6f  =  12.5  Hz. 

As  for  the  current  set  of  measurcniciits,  eight  18  inch  pressure-field  microphones  (G.R.A.S. 
T>"pc  46DD  capsules  with  G.R.A.S.  Type  26CB  preamplifiers)  were  placed  at  grazing  incidence, 
as  shown  in  figure  1.  These  microphones  arc  capable  of  resolving  signals  over  a  broader  range  of 
frequencies  (from  6.5  Hz  up  to  140  kHz  with  ±2  dB  error)  with  a  dynamic  range  from  46  dB  to 
179  dB.  lEPE  power  was  supplied  by  a  Dytran  Instruments  Inc.  4121  power  supply,  which  were 
then  low  pass  hltcrcd  using  a  Khron-Hitc  corp  hltcr  (model  FLX-3007).  These  low  pass  hltcrs 
employ  a  4^^  order  Butterworth  hltcr  and  were  custom  designed  to  have  a  cut  off  frequency  of 
140  kHz.  Microphone  signals  were  digitized  simultaneously  at  500  kS/s  per  channel  using  tvro  NI 
PXI-6122  boards  (four  single  ended  channels  per  board  with  16  bit  resolution  over  ±5v)  for  13.1 
seconds  uninterrupted.  This  largely  over-sampled  rate  ensures  better  resolution  in  the  rise  time  of 
the  shock-structures,  which  is  a  crucial  point  in  this  study.  One  sided  power  spectral  densities  were 
then  generated  from  400  blocks  of  data,  each  comprising  2^"^  samples  per  block  and  a  frequency 
resolution  of  Sf  =  30.52  Hz. 

The  1  8  inch  microphone  measurements  were  conducted  using  four  different  rays  centered  on 
x/Dj  =  17.5,  as  is  shown  in  figure  3.  For  high  Mach  number  supersonic  jets,  the  region  in  the  flow 
between  the  collapse  of  the  potential  core  and  the  supersonic  core  length  arc  known  to  be  the  location 
where  turbulent  mixing  is  not  only  most  intense,  but  also  occurs  at  supersonic  velocities  (15,  16]. 
Identifying  the  correct  propagation  path  is  not  necessarily  trivial  as  the  jet  does  not  behave  like 
a  single  compact  source.  This  is  especially  problematic  close  to  the  hydrodynamic  region  where 
pressure  waves  comprise  a  superposition  of  both  acoustic  and  evanescent  signatures  [17|.  Choosing 
the  wrong  path  artificially  triggers  the  apparition  of  nonlinear  distortions  in  the  signal  that  one  might 
interpret  as  a  sign  of  wave  steepening;  this  problem  has  been  addressed  in  past  studies  |10,  18]. 
Here  we  attempted  to  place  the  microphones  as  close  to  the  jet  as  possible  without  them  being 
endangered,  and  without  contaminating  the  signals  with  evanescent  signatures.  The  growth  of  the 
jet  shear  layer  was  approximated  to  be  O.lx  based  on  the  measurements  of  Tinncy  et  al.  (2008)  |19|. 
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Tabic  2:  Locations  of  1/8  inch  naicrophoncs. 


e 

offset  from  source 

Ap 

65° 

13.0Dj 

lODj 

52° 

iS.ODj 

lODj 

45? 

13.0Dj 

bDj 

455 

35.5Dj 

lODj 

35° 

43.5Dj 

lOD^ 

A  more  conservative  estimate  of  0.15a:  was  employed  as  is  shown  in  figure  3.  Table  2  identifies  the 
location  of  the  1/8  inch  microphone  rays  and  is  listed  in  terms  of  the  ray  angle  (measured  relative 
to  the  jet  centerline),  the  starting  distance  of  each  ray  relative  to  the  post  potential  core  region,  as 
well  as  the  separation  distance  between  subsequent  microphones  on  each  ray. 


Fig.  3:  Location  of  1  8  inch  microphone  rays  and  coordinate  system.  Contour  of  OASPL  (dB, 
ref:  20/iPa^  Hz)  from  14  inch  (x,  r)  grid  measurements  of  Baars  &  Tiiiney  (2014). 

In  an  effort  to  remove  any  possible  contaminations  produced  by  low  frequency  reflections  in 
the  chamber,  a  second  order  high-pass  filter  was  applied  to  the  data  with  a  cut-off  frequency  of 
300  Hz.  Prior  to  this,  corrections  based  on  manufacturer  supplied  data,  was  utilized  to  correct 
for  the  frequency  dependent  error  in  each  transducer.  As  is  shown  in  figure  4a,  the  effect  of  this 
frequency  correction  on  the  spectral  shape  is  subtle,  but  was  used  nonetheless.  Additional  studies 
were  pcrfornuxl  with  and  without  the  microphone’s  protective  caps,  as  well  as  without  the  low- 
pass  antialiasing  filter.  No  discernible  trends  were  observed.  And  so,  both  the  protective  caps  and 
low-pass  filters  were  utilized  for  the  remainder  of  the  studies. 

A  sample  illustration  of  the  final  spectra  is  shown  in  figure  4b,  with  comparison  to  the  FSS  and 
LSS  spectra  of  Tam  et  ai  (1996)  [20].  The  two  microphones  arc  selected  from  the  35°  and  65°  rays 
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Fig.  4:  (a)  One  sided  power  spectral  density  (PSD)  of  a  sample  microphone  signal  at  45°  with  and 
without  corrections  using  manufacturer  supplied  data.(b)  PSD  at  0  =  35°  and  65°  being  compared 
to  the  LSS  and  FSS  of  Tam  et  al.  (1996). 


and  are  loeated  at  85.3Dj  and  83.0Dj  from  the  jet  source.  The  spectra  corresponding  to  the  65° 
observer  is  located  outside  of  the  Mach  cone,  and  so  it  is  characterized  by  a  broader  profile  and  a 
smooth  roll-off,  whereas  the  shallow  angle  obser\^er  at  35°  comprises  a  noticeable  characteristic  peak 
with  relatively  linear  decay.  The  good  agreement  between  the  experimental  data  and  the  similarity 
spectra  is  encouraging  and  demonstrates  the  universality  of  the  spectra  produced  by  high  Mach 
number  jet  flows. 


III.  Sound  pressure  decay  in  the  jet  far-field 
One  of  the  underlying  challenges  in  implementing  the  model  developed  by  Baars  (2013)  [14), 
which  predicts  shock  formation  distances  and  the  degree  of  cumulative  nonlinear  distortion  in  the 
acoustic  waveform  from  high  speed  jets,  is  that  it  requires  one  to  know  the  size  and  amplitude 
of  the  source.  This  is  a  difficult  proposition,  even  with  the  advent  of  more  sophisticated  numer¬ 
ical  simulations  at  our  disposal;  consider  the  different  source  definitions  floating  about  the  open 
literature  [21  23|.  The  experimental  handicap  is  pretty  straightforward:  simply  placing  a  pressure 
sensing  clement  closer  to  a  jet  flow  manifests  a  flow-biased  description  of  the  source  properties,  given 
the  blending  of  both  acoustic  and  hydrodynamic  waves  within  the  jet’s  hydrodynamic  periphery. 
Two  established  methods,  one  described  by  Tinncy  &  Jordan  (2008)  [17],  the  other  by  Grissi  & 
Camussi  (2012)  [24|,  offer  some  relief  in  separating  these  two  components  of  the  pressure  field,  but 
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arc  limited  to  flows  comprising  sound  producing  events  with  eonveetive  speeds  that  arc,  on  average, 
subsonic.  And  so,  an  adequate  means  by  which  the  source  mechanisms,  responsible  for  generating 
Maeh  waves,  can  be  isolated  and  studied  remains  to  be  found. 

A  resurgence  in  recent  years  with  the  use  of  wave  packets  (Morris  (2009)  [25|,  Papamosehou 
(2011)  [26),  Kuo  et  al  (2012)  |27|  and  Jordan  &  Colonius  (2013)  |28|)  may  offer  some  relief,  and 
at  a  reduced  cost.  Modeling  the  evolution  of  supersonically  eonveeting  large-scale  instabilities  by 
wave  packets  provides  a  theoretical  framework  for  studying  the  noise  radiated  by  supersonic  jets. 
In  figure  5  (taken  from  Appendix  A  of  Baars  |14)),  the  standard  deviation  (dp)  of  the  radiating 
component  of  the  pressure  signal  along  the  Maeh  wave  radiation  angle,  and  emanating  from  the 
source  region,  is  shown  to  encompass  three  distinct  zones:  one  that  abides  by  cylindrical  decay 
(cTp  cx  /O'"®  ^),  another  that  abides  by  spherical  decay  (cp  cx  and  a  third  intermediate  zone 

that  transitions  the  two.  Certainly  the  location  where  this  transition  occurs  is  quite  important  in 
determining  source  amplitude.  Baars  (2013)  [14|  shows  how  the  answers  to  this  can  be  expressed 
as  a  ratio  of  the  wave  packet  model  scale  (L)  relative  to  the  jet  diameter,  sueli  that  A  =  L/Dj\  the 
effeet  of  A  on  the  three  decay  regions  is  shown  in  figure  5  for  a  source  located  on  the  lip  line  of 
the  nozzle.  Thus,  to  determine  the  proper  values  associated  with  both  the  wave  packet  scale  and 
its  location  {r^p)-,  a  seiLsitivity  study  was  coiiduetcd.  For  this,  the  spatial  limits  of  the  Gaussian 
envelop  for  all  eases  was  confined  to  three  standard  of  deviations  while  the  wave  packet  amplitude 
was  adjusted  using  pressure  levels  measured  halfway  along  the  array.  Wave  packet  parameters  r.jjp 
and  A  that  were  finally  chosen  were  based  on  a  least  error  sum  between  aU  available  microphones 
and  the  wave  packet.  The  error  function  expression  for  this  is  defined  as  follows, 


(1) 


where  Wi  corresponds  to  the  wave  packet  pressure  \'aluc  located  at  the  nii  microphone  and  i  increases 
with  increasing  distance  from  the  source.  Purthermorc,  measurements  closer  to  the  jet  were  weighted 
heavier  (using  N—i)  than  measurements  further  away,  given  the  larger  ehauges  in  pressure  amplitude 
that  oeeur  closer  to  the  source  region.  The  results  of  this  sensitivity  study  arc  shown  in  figure  6 
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where  the  location  of  the  source  is  expressed  by  its  position  along  the  propagation  path  (po). 
Additional  lines  have  been  added  which  identify  the  final  wave  packet  parameters  {A  —  8.75  and 
Pq  =  (\/2f2)Dj  )  that  were  selected  for  the  remainder  of  the  analysis,  as  well  as  the  location  of  the  jet 
shear  layer.  Thus,  the  large-scale  instability  waves  arc  assumed  to  initiate  their  development  at  the 
exit  plane  of  the  nozzle  and  grow  along  the  nozzle  lip  line  to  their  peak  level  at  x/Dj  =  17.5  before 
decaying  to  negligible  levels  far  downstream.  Under  these  assumptions,  2L  =  17.5.  In  figure  7, 
the  measured  pressure  decay  (using  both  the  1  '4  inch  and  18  inch  microphones  placed  along  the 
45°  path)  compares  favorably  well  to  the  wave  packet  model.  Deviations  from  the  cylindrical  and 
spherical  decay  laws  are  found  to  be  within  2%  of  the  wave  packet  model  at  p/Dj  =  12.3  and  33.4, 
respectively.  The  findings  are  a  testament  to  both  the  scale  of  the  wave  packet  chosen,  as  well  as 
the  propagation  path  along  which  the  microphones  arc  placed. 


Fig.  5:  (a)  schematic  of  the  wave  packet  scale  versus  jet  scale,  (b)  Regions  where  the  wave  packet 
model  obeys  cylindrical  and  spherical  spreading  along  the  peak  noise  path.  Contour  lines  define 

error  thresholds  of  2%,  3%  and  4%. 


Figures  8a  and  8b  present  the  decay  and  inverse  decay,  respectively,  measured  using  the  1/8  inch 
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Fig.  6:  Integrated  errors  between  the  wave  packet  model  and  laboratory  measurements  for  various 
wave  packet  length  scales  and  source  locations  along  the  propagation  path.  Contour  levels  have 
been  normalized  so  that  they  range  between  0  and  1  with  increments  of  0.125. 


Fig.  7:  Spreading  trends  along  the  peak  noise  emission  path  obtained  from  the  wave  packet  model 
using  A  =  8.75.  Both  experimental  datasets  are  shown  with  □  corresponding  to  the  1/4  inch 

microphones. 

microphones  along  all  four  rays.  The  agreement  between  the  45°  and  52°  rays  is  quite  remarkable 
and  confirms  how  linear  decay  holds  well  o\'er  a  range  of  angles  centered  on  the  Mach  wave  angle. 
Thus  one  should  expect  some  flexibility  in  both  the  convective  speed  of  the  relevant  sound  producing 
events,  as  well  as  the  acoustic  region  where  the  propagation  path  represents  that  of  a  wave-guide  so 
to  speak;  spherical  decay  in  the  acoustic  pressure  is  not  restricted  to  a  single  angle.  Albeit,  the  45° 
ray  is  the  only  arrangement  where  a  higher  spatial  resolution  permits  a  verification  that  this  trend 
holds  closer  to  the  jet.  Not  so  surprisingly  is  the  fact  that  the  rate  of  decay  at  steeper  (65°)  and 
shallower  (35°)  angles  is  much  greater;  the  65°  array  being  much  more  extreme.  Glancing  at  the 
OASPL  contours  in  figure  3  justifies  why  this  is  expected. 

A  further  observation  is  made  with  the  trends  presented  in  figures  7  and  8.  Foremost,  the 
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(a) 


Pi^i 

(b) 


Fig.  8;  (a)  cTp  decay  and  (b)  cTp  ^  x  10^  from  1/8  inch  microphones. 


m 


0 


Fig.  9:  OASPL  (dB,  ref:  20/iPa^  Hz)  from  (a;,  r)  plane  grid  incasureincnts  of  Baars  &  Tinney 
(2014)  projected  to  an  are  array  at  p/Dj  =  95  using  spherical  decay.  Assuming  source  region  at 
(a)  x/Dj  =  0  and  (b)  x/Dj  =  17.5.  The  values  shown  are  less  lOOdB. 


remarkable  agreement  between  the  measurements  acquired  with  the  14  inch  and  1  8  inch  iiii- 
erophones  demonstrates  how  transducer  size  and  atmospheric  absorption  have  little  effect  on  the 
pressure  variance  (expressed  here  as  OASPL).  And  so,  the  topography  of  the  OASPL  in  figure  3  us¬ 
ing  1/4  inch  niierophoiies  is  considered  an  accurate  mapping  of  the  pressure  intensity  in  the  far- field 
of  this  Maeh  3  jet  flow.  Likewise,  because  the  sound  pressure  is  shown  to  decay  spherically  after 
approximately  35  diameters,  any  prediction  of  the  far-field  pressure  using  measurements  outside  of 
this  location,  and  assuming  spherical  decay,  should  provide  an  accurate  mapping  of  the  OASPL 
in  the  far-field,  so  long  as  the  correct  propagation  path  is  employed.  It  is  important  that  we  em¬ 
phasize  how  the  location  in  which  the  acoustic  pressure  abides  by  spherical  decay  is  based  on  the 
operating  conditions  of  the  nozzle  and  is  expected  to  scale  with  acoustic  Maeh  number.  Most  often 
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the  propagation  path  is  chosen  to  emanate  from  the  jet  exit  plane.  Baars  &  Tinney  (2014)  [6] 
shows  how  various  features  in  the  far-field  of  supersonic  jets  collapse  exceptionally  well  when  the 
prominent  source  region,  located  after  the  collapse  of  the  potential  core,  is  chosen.  A  demonstration 
of  this  is  shown  in  figure  9  where  the  OASPL  from  the  entire  (x,r)  grid  measurements  of  Baars 
&  Tinney  (2014)  [6]  (120  measurement  points  in  total)  arc  projected  to  an  arc  array  situated  at 
p/Dj  =  95  using  spherical  decay.  In  figure  9a  the  source  is  eissumcd  to  reside  at  the  nozzle  exit 
plane  at  x/Dj  =  0,  while  in  figure  9b  the  source  is  moved  downstream  to  x/Dj  =  17.5.  The  findings 
arc  remarkable  and  expose  the  kinds  of  gross  errors  that  will  occur  when  one  chooses  to  propagate 
jet  noise  data  along  rays  emanating  from  the  jet  exit  plane;  nearly  identical  observations  have  been 
made  in  the  rocket  noise  community  [15,  29|  .  The  only  discrepancies  in  figure  9b  arc  with  the  closer 
measurement  points  (identified  as  filled  and  open  circles)  that  were  located  at  p/Dj  <  35,  and  is 
anticipated  since  these  waveforms  will  initially  undergo  cylindrical  decay  (as  shown  in  figure  7), 
which  was  not  corrected  for  in  this  simple  exercise.  Thus,  both  fine-scale  similarity  and  large-scale 
similarity  like  regions  in  the  far-field  can  be  predicted  from  measurements  closer  to  the  jet  itself,  so 
long  as  the  correct  propagation  path  is  employed. 

IV.  Analysis  of  acoustic  waveforms  along  the  peak  emission  path 
Baars  &  Tinney  (2014)  |6|  flenionstratcd  how  the  statistical  properties  of  the  acoustic  waveforms 
are  relatively  unchanged  along  rays  emanating  from  the  prominent  somce  regions  of  supersonic  jets; 
similar  results  were  shown  here  in  section  III.  Furthermore,  only  weak  waveform  distortions  were 
measured  using  1/4  inch  microphones  placed  in  the  mid  to  far-field  regions  of  the  jet  and  under 
higher  levels  of  humidity.  And  so,  the  motivation  of  this  part  of  our  study  is  twofold.  The  first  is 
to  demonstrate  the  effect  of  sensor  size  and  atmospheric  absorption  on  capturing  these  distortions, 
while  the  second  is  to  address  the  postulations  of  Baars  k,  Tinney  (2014)  |6]  through  a  thorough 
analysis  of  the  data  acquired  with  the  higher  fidelity  instruments. 

The  statistical  metrics  that  are  employed  comprise  the  quadrature  spectral  density,  skewness 
of  the  pressure  time  derivative,  the  number  of  zero  crossings  per  unit  time  [Nze)  and  the  integral 
of  the  negative  component  of  the  quadrature  spectral  density.  These  metrics  have  been  thoroughly 
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described  elsewhere  and  arc  used  to  determine  the  degree  of  nonlinearity  in  the  acoustic  waveform 
at  a  single  point.  When  analyzed  over  several  points  along  the  propagation  path,  the  changes 
to  the  waveform  amount  to  nonlinear itics  cncumberefl  from  eumulativc  distortions.  This  is  true 
so  long  as  the  measured  waveforms  are  from  a  single  source,  or  have  limited  interference  from 
neighboring  wavefronts  produced  by  other  nearby  sources  in  the  flow.  The  latter  will  give  rise  to 
artificial  steepening  and  the  formation  of  rogue  waves  formed  from  coalescence.  Other  metrics  were 
considered  and  included  the  wave  steepening  factor  and  Kurtosis,  though  neither  revealed  anything 
significant  and  so  they  arc  not  shown  here.  A  listing  of  these  metrics  is  provided  as  follows: 

•  Quadrature  Spectral  Density.  The  definition  for  this  metric  is  well  documented  in  the  work 
of  Morfey  &  Howell  (1981)  [30]  and  Howell  &  Morfey  (1987)  [31  j.  This  quantity  highlights  the 
importance  of  non-linear  distortion  across  the  spectrum  relative  to  the  power  sj)cctral  density; 
the  normalized  forms  of  this  indicator  is  as  follows, 


(2) 


Q  _  Qp^p/^p  .  . 

5  Spp/CTI 

where  Qp2p  =  —  ^(2p2(/)p*(/)),  and  P2{f)  and  p(f)  are  the  Fourier  transforms  of  the  signals 
p{t)  and  p^(t),  respectively.  The  denominator  in  Eq.  (3)  comprises  the  ensemble  averaged 
power  spectral  density,  which  is  simply  5pp  =  {p{f)p*(f))^  Howell  &  Morfey  (1987)  |31| 
recognized  early  on  the  limitation  in  using  these  metrics  at  a  single  point  and  that  the  problem 
[of  predicting  distortion]  remains  unsolved  without  estimating  Qp2p  as  a  function  of  distance. 
This  is  especially  problematic  with  experimental  data,  where  transducer  errors  inherently 
prevent  the  signal  from  being  truly  Gaussian;  thus  it  is  difficult  to  envision  this  term  ever 
being  zero  in  a  laboratory  environment.  Figure  8a  illustrates  how  the  pressure  standard 
deviation  along  the  45°  ray  abides  by  realistic  decay  laws,  and  so  an  application  of  both  Eqs. 
(2)  and  (3)  to  this  data  should  provide  meaningful  results. 
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•  Skewness  of  the  pressure  derivative.  The  amount  of  wave  steepening,  and  thus  the  degree 
of  nonlinearity  in  the  pressure  signal,  ean  be  quantified  by  eomputing  the  probability  density 
funetioii  (PDF)  and  skewness  faetor  of  the  pressure  time  derivative  p(t).  The  PDF  is  denoted 
as  D  (p)  while  its  asymmetry  is  defined  by  the  third  eentral  moment. 


(4) 


The  dimensionless  asymmetry  is  thus  the  well  known  skewness  factor,  and  is  defined  as 


where  dp  is  the  variance  ofp(i).  Numerous  applications  using  this  metric  to  quantify  nonlin- 
earities  in  real  waveforms  ean  be  found  in  the  open  literature  with  the  skewness  of  the  pressure 
time  derivative  providing  the  most  sensitivity  [32|. 

•  Number  of  zero  crossings  per  unit  time.  Since  higher  amplitude  waves  move  faster 
along  the  propagation  path  and  merge  or  coalesce  with  weaker  waveforms,  the  number  of  zero 
crossings  will  decrease  with  increasing  distance  from  the  source.  This  does  not  necessarily 
mean  that  the  waveform  is  steepening  along  the  propagation  path  and  so  the  link  between 
Nzr  and  cumulative  nonlinear  distortion  effects  warrants  caution. 

•  Integral  of  negative  Q,  If  cumulative  distortions  are  occurring  along  the  propagation  path, 
then  one  would  anticipate  a  flux  of  energy  in  and  out  of  the  mid  frequencies.  This  ean  be 
quantified  by  integrating  the  negative  part  of  the  quadrature  spectral  density  at  each  location 
on  the  ray  as  follows, 


Q<0 
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A.  Spectral  indicators 

Profiles  of  Q  and  Q/S  are  shown  in  figures  10a  and  10b  using  equations  (2)  and  (3),  respectively, 
along  the  four  different  rays  (a  1%  bandwidth  moving  filtered  has  been  employed).  As  expected,  both 
Q  and  Q/S  profiles  reflect  a  loss  in  energy  at  mid  frequencies,  with  Q/S  demonstrating  an  increase 
in  energy  at  high  frequencies.  Closer  inspection  reveals  two  distinct  trends.  The  first  comprises  grey 
lines  with  arrows  indicating  the  direction  of  increasing  distance  from  the  source  and  with  noticeable 
changes  in  both  energy  and  peak  frequency.  Increases  in  mid-frequcncy  energy  suggest  waveform 
coalescence  while  increases  in  high  frequency  energy  suggests  waveform  steepening.  The  other  set, 
highlighted  by  solid  dark  lines,  encompasses  regions  that  reside  beyond  approximately  40Dj.  These 
black  lines  collapse  reasonably  well  and  suggest  that  the  acoustic  waveforms  are  no  longer  coalescing 
or  distorting  along  the  propagation  path;  this  agrees  well  with  the  measurements  of  Baars  Tiimey 
(2014)  [6|  using  the  14  inch  microphones.  Along  the  65°  ray,  the  trends  are  difficult  to  decipher. 
Since  this  ray  is  outside  the  Mach  cone,  the  lack  of  any  discernible  trend  is  expected.  Albeit,  the 
near  perfect  collapse  of  the  OASPL  shown  in  figure  9b  is  a  testament  to  the  choice  of  propagation 
path,  even  for  the  sideline  and  steep  angle  observer.  This  is  consistent  with  the  measurements  of 
Baars  Sz  Tinney  (2014)  [6]  where  it  has  been  shown  that  shock  structures  are  concentrated  mainly 
along  the  Mach  wave  angle. 

B.  Scalar  indicators 

Figure  11  compares  the  number  of  zero  crossings  per  second,  the  skewness  of  the  pressure 
derh^ative  and  the  integral  of  Qnrg  along  the  45°  path.  The  data  acquired  with  the  1  8  inch 
microphones  are  analyzed  using  both  the  original  sample  rate  of  500  kHz  as  well  as  a  reduced  rate 
of  100  kHz,  which  coincides  with  the  sample  rate  used  by  the  1  /4  inch  microphones.  The  findings 
from  the  latter  are  used  to  compare  the  1/8  inch  microphone  measurements  with  the  data  acquired 
using  the  1  '4  inch  microphones  in  an  effort  to  determine  the  effect  of  transducer  size  and  atmospheric 
properties  (humidity  and  temperature)  on  these  higher  order  metrics.  For  each  metric,  an  arbitrary 
approximation  of  the  trend  corresponding  to  the  500  kHz  sampling  of  the  1/8  inch  inicrophone  is 
shown,  along  with  its  spatial  derivative.  These  deri\^tives  reveal  three  distinct  zones.  The  first 
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(a)  (b) 

Fig.  10:  (a)  Q  x  10“^  and  (b)  Q/5,  Morfcy-Howell  indicators  (1%  bandwidth  moving  filter)  for  all 
microphones  (gray)  and  for  the  last  4  microphones  (black). 


manifests  a  fairly  constant  spatial  rate  of  change  while  the  second  reveals  changes  in  amplitude 
along  the  propagation  path.  The  third  zone  identifies  a  region  where  the  metric  is  saturated  and  is 
no  longer  changing. 

Where  the  number  of  zero  crossings  and  the  skewness  of  the  pressure  time  derivative  are  con¬ 
cerned,  the  effect  of  sampling  rate  is  clearly  manifest.  However  the  effects  of  temperature  and 
humidity  on  the  number  of  zero  crossing  and  the  integral  of  Qneg  appear  not  to  be  affected.  While 
corrections  for  atmospheric  absorption  arc  well  established,  they  cannot  be  corrected  in  the  raw  time 
signal.  And  so  the  effects  of  atmospheric  absorption  can  only  be  ascertained  from  raw  measurements 
of  this  kind. 

As  for  the  skewness  of  the  pressure  derivative  7(p),  several  relevant  features  are  seen  in  fig¬ 
ure  lib.  First  of  all,  three  trends  corresponding  to  each  data-set  (same  legend  as  in  figure  11a) 
are  clearly  present.  The  new  data-set  sampled  at  500  kHz  shows  a  quasi-lincar  growth  of  7(p)  that 
suggests  wave  steepening;  zero  skewness  is  reflective  of  a  normal  distribution.  By  80  or  so  jet  di- 
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amctcrs  (/o  =  2  m)  the  steepening  process  appears  to  have  reached  its  niaximuin  at  which  point  the 
trend  remains  either  constant  or  begins  to  decay.  For  the  resampled  data  (1/8  inch  microphones 
resampled  at  100  kHz),  the  change  in  7(p)  follow’s  the  originally  sampled  data  consistently  up  until 
roughly  0.7  meters  from  the  source  region,  after  which  changes  in  the  skewmess  begin  to  reflect  the 
trends  produced  by  the  1  M  inch  microphone  study.  This  agreement  demonstrates  how  distortion 
mechanisms  affect  the  low  frequency  bands  in  the  w'avcform  earlier  on  2is  the  lower  sampled  data 
reflects  that  of  the  higher  sampled  data.  Differences  between  the  resampled  18  inch  microphone 
data  and  the  1/4  inch  microphone  measurements  arc  attributed  to  elevated  levels  of  humidity  during 
the  1/4  inch  study.  A  visual  extrapolation  of  the  trends  toward  the  origin  {p  — >  0)  suggests  that 
the  waveforms  produced  in  the  source  region  of  the  jet  arc  already  skewed,  thus  confirming  the 
simulations  of  Nichols  et  a/.(2013)  [ll|. 

Figure  lie  quantifies  the  net  flux  of  energy  transfer  across  the  quadrature  spectral  density 
shown  in  figure  10a.  By  restricting  the  limits  of  integration  to  only  the  negative  part  of  the  spectra, 
which  is  showm  to  be  concentrated  on  the  mid  frequency  range,  any  noise  inherent  to  the  borders 
of  the  spectrum  arc  avoided.  This  registers  the  flux  of  mid-range  spectral  energy  that  is  lost  to 
w’avcform  coalescence  and  is  therefore  a  relevant  indicator.  Once  again,  a  trend  very  similar  to  7(p) 
appears  with  a  noticeably  linear  growth  at  the  start  of  the  array.  Further  away  the  evolution  starts 
to  decrease  and  eventually  converges  as  the  spatial  rate  of  change  becomes  negligible. 

The  evolution  of  the  spatial  rate-of-cliangc  along  the  array  for  these  three  indicators  arc  shown 
alongside  one  another  in  figure  12  with  each  one  being  normalized  by  its  maximum  amplitude  in 
order  to  compare  trends.  Overall,  these  indicators  comprise  a  similar  shape,  which  w^ould  suggest 
that  both  waveform  stecpciiiiig  and  coalescence  are  occurring  within  2.5  meters  or  so  from  the  source 
region  of  this  laboratory  scale  jet. 

C.  Shock  formation  distance 

In  an  effort  to  link  the  observed  phenomena  to  the  theoretical  estimates  for  shock  formation 
distance,  two  such  formulations  arc  ooiisidercd,  one  for  spherically  propagating  waves  and  the  other 
for  cylindrical  weaves.  Given  the  changes  in  pressure  amplitude  registered  by  our  measurements 
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Fig.  11;  (a)  Number  of  zero-crossings,  (b)  skewness  of  the  pressure  time  derivative  and  (c) 

—Qneg  X  10^  along  the  45°  ray. 
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Fig.  12:  Comparison  of  the  spatial  derivative  of  Nze^  7(7^)  and  Qneg  normalized  by  their  respective 

maximum  value. 


and  confirmed  by  wave  packet  modeling,  the  real  solution  will  be  a  blended  cstiiiiatc  of  the  tw^o. 
Following  Hamilton  k.  Blackstock  (2008)  [8|,  the  shock  formation  distance  for  an  initially  broadband 
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gaussiaii  wavcfonii  that  propagates  cylindrically  and  spherically  from  its  source  in  a  lossless  fluid 
can  be  estimated  from,  respectively, 

rc  =  ro(l  +x/2ro)^  (7) 

and 

fs  =  rocxp(x/ro),  (8) 

where  x  is  the  plane  wave  shock  formation  distance  defined  by  Gurbatov  &  Rudenko  (p.383  of 
Hamilton  Sz  Blackstock  [8])  as, 

/3(27r/o)<Tp,o’  ^  ^ 

and  ro  is  the  source  radius.  The  coefficient  of  nonlinearity  is  given  by  =  (7  +  l)/2  =  1.2  while 
the  density  of  air  is  determined  from  table  1  and  ideal  gas  properties  to  be  =  1-14  kg/m^. 
Figure  4b  reveals  a  cliaractcristic  frequency  valued  at  fo  =  3.G  kHz  for  the  45°  ray,  which  is  the 
same  used  by  Baars  (2013)  [14].  Both  the  size  and  strength  of  the  source  arc  not  easily  determined, 
given  the  expcriincntal  and  numerical  handicaps  discussed  in  section  III.  Likewise,  expressions  for 
shock  formation  distance  require  that  A:r  ^  1,  that  is,  the  source  radius  must  exceed  the  maximum 
wavelength  of  the  radiated  sound.  Three  likely  locations  for  the  source  surface  are  considered  in  this 
exercise,  though  an  infinite  number  of  solutions  can  be  obtained.  The  first  borders  the  periphery 
of  the  jet  shear  layer  at  p  =  3.5Dj  by  assuming  a  shear  layer  growth  of  O.lx  [19|.  The  second  is 
selected  to  replicate  the  location  used  by  Baars  (2013)  [14]  so  that  ro  is  located  at  p  =  2.5Dj.  The 
third  location  is  assumed  to  be  located  on  the  nozzle  lip-line  at  p  =  (\/2/2)Dj,  where  the  peak  noise 
source  has  been  shown  to  reside  [35].  Each  of  these  locations  is  within  the  cylindrical  decay  region  of 
the  pressure  amplitude,  and  so  the  source  amplitudes  (for  these  three  points)  arc  extrapolated  using 
the  wave  packet  solution  for  L  =  8.75  and  assuming  p  oc  l/>yp.  Shock  formation  distances  arc  then 
estimated  using  Eqs.  (7)  and  (8)  for  cylindrical  and  spherical  decay  along  the  entire  propagation 
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Tabic  3:  Shock  formation  distances  and  effective  Gordberg  numbers. 
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path,  respectively.  Thus,  for  a  starting  waveform  that  abides  by  gaussian  statistics  and  is  allowed  to 
propagate  in  the  absence  of  viscous  absorption  effects,  its  shock  formation  distance  will  be  bounded 
by  Tc  and  Ts . 

Findings  from  this  analysis  arc  shown  in  table  3  alongside  the  original  estimates  from  Baars 
(2013)  |14|  who  estimated  the  source  amplitude  to  be  a  factor  of  three  larger  than  what  is  measured 
here;  Baars  (2013)  ll4|  assumed  the  source  amplitude  abided  by  spherical  propagation  along  its  entire 
path,  which  is  reasonable,  considering  the  information  that  was  available  at  the  time.  Nonetheless, 
the  recipe  outlined  by  Baars  (2013)  [14|  is  employed  here  and  provides  reasonable  estimates  for 
Tr  and  r^.  The  findings  suggest  that,  aside  from  both  the  effects  imposed  by  having  presteepened 
waveforms  at  the  source  as  well  as  viscous  relaxation  along  the  propagation  path,  it  is  unlikely  that 
s'hoeks  will  form  within  the  eonfincs  of  this  laboratory  facility. 

D.  Gol’dberg  number 

Naturally,  the  next  step  is  to  calculate  the  Gol’dbcrg  number,  a  measure  of  the  strength  of 
nonlinear  distortion  relative  to  that  of  dissipation.  When  F  <  1,  attenuation  dominates  and  the 
formation  of  shocks  is  suppressed.  W’hen  F  1  nonlinear  effects  will  be  distinctly  present.  The 
expression  for  plane  waves  is  well  established  and  is  defined  as, 

r  =  '^  (10) 

X 

where  the  absorption  length  scale  (ia)  is  simply  the  inverse  of  the  absorption  coeffieient  (a)  which 
is  evaluated  in  table  3  using  an  expression  from  Bass  et.  al.  (1996)  [36].  W^e  evaluate  a  using  the 
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characteristic  frequency  /o  =  3.6  kHz  and  for  a  relative  humidity  determined  by  the  average  value 
recorded  during  the  45°  2  tests  (%RH=  45).  The  resultant  absorption  for  the  1/8  inch  microphone 
study  is  found  to  be  a  =  0.0026  Np/m,  which  yielded  £a  =  385  m  {ia  —  417  m  for  the  1  '4  inch 
microphone  study  conducted  by  Baars  et  al.  (2012)  [10]  based  on  70%  relative  humidity  and  an 
atmospheric  temperature  of  290  K).  Surprisingly,  the  differences  in  absorption  between  the  current 
measurements  and  the  1/4  inch  study  is  small;  changes  in  temperature  are  balanced  by  changes  in 
humidity  for  this  particular  frequency.  As  was  done  for  the  shock  formation  distance,  an  expression 
for  the  Gol’dberg  number  for  cylindrically  and  spherically  propagating  waves  is  needed.  Hamilton 
(2013)  |37|  recently  showed  how  the  effective  GoLdberg  number  for  diverging  waves  can  be  expressed 
by  a  single  parameter  as  follows, 

”  1  +  7r^/4ro 

A.S  =  rcxp(— 7rx/2ro)  (12) 

where  Eq.  (11)  and  (12)  arc  the  cylindrical  and  spherical  forms  of  the  expression.  Estimates  of  Ac 
and  As  for  the  1  8  inch  and  1  4  inch  studies  are  provided  in  table  3  with  contradictory  conclusions 
being  offered.  That  is,  estimates  that  assume  purely  cylindrical  decay  suggest  cumulative  nonlinear 
distortions  will  be  present  in  the  waveform,  whereas  estimates  assuming  pure  spherical  decay  suggest 
that  no  such  distortions  will  form. 

E.  Direct  observation 

A  comparison  of  the  instantaneous  waveforms  registered  along  the  45°  ray  in  the  inner  and  outer 
zones  arc  shown  in  figures  14  and  15,  respectively.  These  waveforms  are  aligned  by  cross-eorrclating 
the  first  microphone  signal  with  subsequent  signals  on  the  array,  where  ^pp'ir)  =  {p{p,t)p{p\t)). 
The  corrected  time  delay  is  identified  by  the  maximum  correlation,  as  shown  in  figure  13,  which 
exposes  two  important  features.  That  is,  correlations  of  this  magnitude  demonstrate  how  directive 
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the  waveforms  arc  along  this  path,  while  the  decay  in  correlation  level  with  increasing  distance  is 
evidence  of  contamination  from  other  waveforms  produced  by  other  sources  located  elsewhere  in 
the  flow.  And  so,  this  raises  additional  concerns  as  to  whether  cumulative  waveform  distortions 
are  forming  beyond  the  hydrodynamic  near-field  regions  of  laboratory  scale  jets,  aside  from  the 
statistical  analysis  provided  in  Section  IV  A. 


T  •  aoc/Dj  +  35.5 

Fig.  13:  Cross  corrections  between  inicrophones  on  the  45^  ray  in  the  outer  zone  eiicoiiipassing 

p/Dj  =  [35.5  :  10  :  105.5]. 

Turning  our  attention  to  the  raw  pressure  time  series  in  figures  14  and  15,  several  features 
arc  manifest,  which  complement  the  simulations  of  Nichols  et  al.  (2013)  [ll|  and  Anderson  & 
Freund  (2012)  [12].  Albeit,  without  knowing  the  particle  velocity  simultaneously,  the  discussion 
and  interpretation  provided  here  warrants  caution.  In  figures  14  and  15,  several  additional  lines 
have  been  drawn  with  labels  corresponding  to  the  following  kinds  of  events: 

•  Wave  Steepening  (WS).  This  corresponds  to  instances  when  the  wave  registered  at  multiple 
locations  on  the  ray  appears  to  be  the  same  wavefront  and  is  steepening.  This  is  identified  by 
a  single  wave  that  evolves  from  a  more  relaxed  wave  shape  to  a  N-wave  shape  with  increasing 
distance. 

•  Wave  Passing  (WP).  These  waves  correspond  to  instances  where  the  direction  of  propagation 
differs  slightly  from  the  45°  path  along  which  the  microphones  arc  plaeed.  Because  the  width  of 
the  wavefront  is  finite,  it  either  appears  or  disappears  in  the  microphone  signal  at  various  points 
consecutively.  It  is  assumed  that  the  number  of  waves  entering  and  exiting  the  propagation 
path  is  the  same. 

•  Rogue  Wave  (RW).  These  will  form  when  two  smaller  waveforms  travel  at  small  angles  relative 
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to  one  another  and  eventually  merge  to  produce  a  much  larger  and  steepened  wave  whose 
amplitude  is  approximately  that  of  the  original  waves  combined. 


•  Wave  Coalescence  (WC).  This  would  appear  similar  to  a  rogue  wave  where  two  waves,  prop¬ 
agating  along  the  same  path,  merge  and  then  coalesce. 

The  significance  of  having  both  passing  waves  and  rogue  waves  occurring  along  the  propagation 
path  is  that  both  will  artificially  escalate  Nze^  7(p)  and  Qneg  and  mislead  one  into  believing  that 
cumulative  distortions  are  occurring  along  the  propagation  path.  For  the  sample  time  series  shov;n 
in  figures  14  and  15,  instances  of  passing  waves  is  significant.  Rogue  waves  only  appear  in  these 
time  series  at  further  distances  in  figure  15  where  N-wave  like  structures  are  also  significant. 
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Fig.  14:  Evolution  of  the  raw  pressure  time  series  (shown  as  p{p,t)  j ap{p))  along  the  45°  ray  in  the 
inner  zone  encompassing  p/Dj  =  [13  :  5  :  48). 


V.  Conclusions 

The  acoustic  waveforms  produced  by  an  unheated  Mach  3  jet  are  investigated  using  1  4  inch 
microphones  covering  a  2-D  grid  in  the  (a:,r)  plane  as  well  as  1/8  inch  microphones  placed  along  rays 
emanating  from  the  source  regions  of  the  flow  at  x/Dj  —  17.5.  Various  metrics  including  pressure 
amplitude  decay,  the  normalized  quadrature  spectral  density,  skewness  of  the  pressure  derivative, 
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Fig.  15:  Evolution  of  the  raw  pressure  time  series  (shown  as  p(p,  i)/<Tp(p))  along  the  45^  ray  in  the 
outer  zone  encompassing  p/Pj  =  [35.5  :  10  :  105.5]. 


the  number  of  zero  crossings  and  the  integral  of  the  negative  part  of  the  quadrature  spectral  density 
were  used  to  study  the  behavior  of  these  waveforms  along  the  propagation  path.  Where  the  pressure 
amplitude  is  eoneerned,  the  location  where  its  rate  of  decay  first  abides  by  cylindrical  spreading, 
followed  by  spherical  spreading,  was  identified  and  was  confirmed  using  a  simple  wave  packet  model. 
The  OASPL  appeared  not  to  be  affected  by  sensor  size  and  so  the  linear  prediction  of  its  amplitude 
at  95  jet  diameters  was  performed  and  shown  to  collapse  remarkably  well  for  all  azimuthal  angles 
when  propagated  along  rays  emanating  from  the  post  potential  core  regions  of  the  jet. 

Both  spectral  and  scalar  metrics  used  to  eonipute  nonlinear  waveform  distortions  revealed  three 
different  areas  in  the  jet  near-field.  The  first  comprised  a  zone  where  waveform  distortions  appear 
to  be  increasing  at  a  constant  rate  that  initiate  in  close  proximity  to  the  jet  (within  one  meter  from 
the  prominent  source  region).  A  second  zone  showed  a  decrease  in  nonlinear  distortion  effects  while 
a  third  zone  demonstrated  how  the  waveforms  were  no  longer  steepen  and  had  either  reached  their 
maximum  or  were  beginning  to  undergo  a  state  of  relaxation.  A  visual  extrapolation  of  the  nonlinear 
metrics  towards  the  source  region  suggests  that  the  acoustic  waveforms  are  initially  distorted  at  the 
source.  Consequently,  a  direct  observation  of  the  raw  acoustic  time  scries  along  the  propagation 
path  revealed  both  passing  waves  and  rogue  waves  in  the  pressure  waveform,  both  of  which  would 
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ultimately  escalate  the  metrics  for  quantifying  cumulative  nonlinear  distortions. 
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Appendix  B:  2014  Texas  Fluid  Dynamics  IMeeting  Presentation 
Nonlinear  Acoustic  Distortion  in  Supersonic  Lab-Scaled  Jet  Vicinity 
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Texas  Fluid  Dynamics  Meeting,  2014 

Nonlinear  acoustic  distortion  in  supersonic  lab-scaled  jet  vicinity. 
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The  spatial  evolution  of  acoustic  waveforms  produced  by  a  Mach  3  jet  are  investigated  in  a  fully  anechoic 
chamber  using  both  1/4  inch  and  1/8  inch  pressure  field  microphones  located  along  rays  emanating  from 
the  post  potential  core  where  the  peak  sound  emission  is  found  to  occur.  Predictions  of  the  OASPL  along 
an  arc  located  at  95  jet  diameters  using  120  planar  grid  measurements  are  shown  to  collapse  remarkably 
well  when  the  arc  array  is  centered  on  x/Dj  =  17.5  (post  potential  core),  as  opposed  to  the  jet  exit  plan. 
Statistical  metrics  such  as  the  quadrature  spectral  density,  number  of  zero  crossings,  the  skewness  of  the 
pressure  time  derivative  and  the  integral  of  the  negative  part  of  the  quadrature  spectral  density  are  used. 
The  results  demonstrate  how  distorted  waves  form  through  coalescence  and  wave  steepening  and  occur 
within  a  few  meters  from  the  source  in  these  laboratory  scale  jets.  The  remarkable  agreement  between  a 
wavepacket  modeling  of  the  jet  near-field  and  the  experiment  allows  identifying  the  transition  zone  from 
cylindrical  to  spherical  pressure  decay  law,  and  demonstrates  how  the  latter  of  the  two  occurs  within  35 
Dj  along  the  propagation  path. 
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